Ab initio nuclear structure simulations: the speculative 14 F nucleus 
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We present results from ab initio No-Core Full Configuration simulations of the exotic proton-rich 
nucleus 14 F whose first experimental observation is expected soon. The calculations with JISP16 
NN interaction are performed up to the iV max = 8 basis space. The binding energy is evaluated 
using an extrapolation technique. This technique is generalized to excitation energies, verified in 
calculations of 6 Li and applied to 14 B, the 14 F mirror, for which some data are available. 



PACS numbers: 21.60.De, 27.20.+n, 21.10.Dr 

Exotic nuclei at the nuclcon drip lines and beyond 
constitute a forefront research area in nuclear physics. 
The physics drivers include: (1) to discover how shell 
structures evolve into extreme isospin regions; and (2) 
to extend our knowledge of the strong interactions, es- 
pecially elusive three- nucleon forces (3NF), under these 
conditions. In order to help pave a path towards these 
goals, we present baseline ab initio results for selected 
unstable A=13 and A=14 nuclei. 

We focus especially on 14 F, with isospin T=2, that is 
expected to lie beyond the proton drip-line and therefore 
unstable. This proton-rich nucleus will strain the conver- 
gence properties of the ab initio methods we adopt here, 
and also push the limits of state-of-the-art experimental 
facilities. Indeed, the first experimental results regard- 
ing this four proton excess isotope will be available soon 
from Cyclotron Institute at Texas A&M University 

We perform the first ab initio study of 14 F. We use 
the No-Core Shell Model (NCSM) @, [| which employs 
a many-body harmonic oscillator basis which treats all 
nucleons as spectroscopically active. The basis space in- 
cludes all many-body states with excitation quanta less 
then or equal to N max that makes it possible to com- 
pletely remove spurious center-of-mass excitations. Wc 
used the code MFDn 0-Q and the realistic NN inter- 
action JISP16 1 @ which is known to provide a good de- 
scription of p shell nuclei 0, H| without an additional 
3NF. The largest calculations were performed in the 
Njoax = 8 basis space, which for this nucleus contains 
1,990,061,078 basis states with total magnetic projection 
M = and natural parity (negative). The determina- 
tion of the lowest ten to fifteen eigenstates of the sparse 
Hamiltonian matrix, for each oscillator parameter Ml, re- 
quires 2 to 3 hours on 7,626 quad-core compute nodes at 
the Jaguar supercomputer at ORNL. 

We show our complete results for the 14 F ground state 
energy in Fig. [TJ The solid curves are our NCSM results 
with the bare interaction JISP16. These results are strict 
upper bounds for the ground state energy, and converge 



1 A Fortran code for the JISP16 interaction matrix elements is 
available at http://nuclcar.physics.iastatc.edu 



monotonically with A max to the infinite basis space re- 
sults. The dashed curves in Fig. Q] are obtained in more 
conventional NCSM calculations with effective NN inter- 
actions derived from the initial bare interaction JISP16 
by the Lee-Suzuki-Okamoto (LSO) rcnormalization pro- 
cedure The renormalization procedure is truncated 
at the two-body cluster level — i. e. induced three-body, 
four-body, etc., contributions are neglected; hence we re- 
fer to these calculations as LSO(2) renormalized. Note 
that these results differ slightly from preliminary approx- 
imate results presented at recent conferences [10|, HI! • 

By comparing the bare and LSO(2) renormalized 
JISP16 results in Fig. Q] we observe that the tendency 
of the LSO (2) renormalized calculations is misleading. 
For increasing basis spaces from N max = to 6, the min- 
imum of the ftf2-dependent curves increases, suggesting 
an approach from below to the infinite basis space result. 
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FIG. 1: (Color online) Results for the ground state energy 
of 14 F with bare (solid) and LSO (2) renormalized (dashed) 
JISP16 as a function of the oscillator parameter Ml. The 
shaded area demonstrates the global extrapolation A for the 
binding energy and its uncertainty; the extrapolation B at 
fixed h~Q is given by stars. The most reliable Ml value for 
this extrapolation method is at Ml = 25 MeV for 14 F with its 
uncertainty indicated by the error bar. 



2 



At iV max = 6, LS0(2) renormalizcd JISP16 produces a 
nearly flat region at approximately the same energy as 
the minimum obtained with the bare JISP16 interaction. 
On the other hand, the bare interaction provides a vari- 
ational upper bound for the ground state energy, which 
decreases with increasing iV max . 

Other light nuclei ( 6 He, 6 Li, 8 Be, 12 C, 16 0, . . . ) show 
a qualitatively similar behavior: the LSO renormalized 
interactions produce results which are neither an upper 
bound nor a lower bound, and the approach to the in- 
finite basis space is non-monotonic. Hence the conver- 
gence pattern of the LSO renormalized results is diffi- 
cult to assess. Furthermore, with the patterns displayed 
in Fig. Q] for JISP16, the minima of the ^-dependent 
ground state energy curves for both the bare and the 
LSO(2) renormalized interaction may be expected to co- 
incide for N max > 8 as in some other nuclei. For these 
reasons, we did not perform expensive N max = 8 LSO(2) 
renormalized JISP16 calculations for 14 F. 

Recently we introduced the ab initio No-Core Full Con- 
figuration (NCFC) approach [H, lli^ . by extrapolating 
NCSM results with the bare interaction in successive ba- 
sis spaces to the infinite basis space limit. This makes 
it possible to obtain basis space independent results for 
binding energies and to evaluate their numerical uncer- 
tainties. We use two extrapolation methods: a global 
extrapolation based on the calculations in four succes- 
sive basis spaces and five Ml values in a 10 MeV interval 
(extrapolation A) ; and extrapolation B based on the cal- 
culations at various fixed Ml values in three successive 
basis spaces and defining the most reliable Ml value for 
the extrapolation. These extrapolations provide consis- 
tent results Q. Combining both extrapolation methods 
suggests a binding energy of 72 ± 4 MeV for 14 F which 
agrees well with AME03 nuclear binding energy extrap- 
olations [111, see Tablefl] Ironically, out of all our NCSM 
calculations, both with the bare and the LSO(2) renor- 
malized interaction, the minimum of the LSO (2) calcu- 
lations at V max = appears to be closest to the infinite 
basis space result. 

To check the accuracy of our approach, we performed 
similar calculations for the mirror nucleus 14 B with a 
known binding energy of 85.423 MeV [l3|. This value 
agrees with our NCFC result of 86 ± 4 MeV. We also 
performed NCFC calculations of the neighboring nucleus 
13 O using basis spaces up to N max = 8. The calculated 
binding energy of 77 ± 3 MeV also agrees with the exper- 



TABLE I: NCFC results obtained with JISP16 for the ground 
state energies (in MeV) of 13 0, 14 B and 14 F. Experimental 
data are taken from Ref. [l3l |. 



Nucleus 


Extrap. A 


Extrap. B 


Experiment 




-75.7(2.2) 


-77.6(3.0) 


-75.56(0.01) 


14 B 


-84.4(3.2) 


-86.6(3.8) 


-85.42(0.02) 


14 F 


-70.9(3.6) 


-73.1(3.7) 


-73.3(0.4)° 



a AME03 extrapolation fli . 



imental value of 75.556 MeV [H. 

Wc note that a good description of both 14 F and 13 in 
the same approach is important in order to have a consis- 
tent description of the 13 0+p reaction that produces 14 F. 
In this way, any experimentally observed resonances can 
be directly compared with the difference of our results for 
the 14 F and 13 energies. In this respect it is interesting 
to note that although the total energies of the extrapo- 
lations A and B differ by about 2 MeV, the differences 
between the ground state energies of these three nuclei 
are almost independent of the extrapolation method: for 
14 F and 13 the predicted difference is 4.6 MeV, and for 
14 F and 14 B it is 13.5 MeV. (The numerical uncertainty 
in these differences is unclear, but expected to be signifi- 
cantly smaller than the uncertainty in the total energies.) 

We also calculated the 14 F excitation spectrum in an- 
ticipation of the experimental results. It is unclear how 
to extrapolate excitation energies obtained in finite basis 
spaces, but we can extrapolate the total energies of ex- 
cited states using the same methods as discussed above 
for the ground state energy. For the lowest state in each 
J 71 channel the convergence pattern should be similar 
to that of the ground state; for excited states with the 
same quantum numbers we simply assume the same con- 
vergence pattern. We perform independent separate ex- 
trapolation fits for all states. The differences between the 
extrapolated total energies and the ground state energy 
is our prediction for the excitation energies. 

This approach to extrapolating the total eigenenergies 
is supported by applying it to 6 Li, see Table [Til We have 
results for 6 Li in basis spaces up to V max = 16 where a 
good convergence is achieved and hence the extrapolation 
uncertainties are small. These results are compared in 
Table [TTJ with the extrapolations based on calculations 
in basis spaces up to N max = 8, i. e. in the same basis 
spaces used for the 14 F and 14 B studies. 

We see that the excitation energies based on N max = 8 
and smaller basis space results are consistent with the 
results obtained in larger spaces. The level ordering is 
the same and the difference between the V max = 8 and 
Nmax = 16 results is generally much smaller than the esti- 
mated uncertainties in the total energies of the N max = 8 
extrapolations. This suggests that the numerical uncer- 
tainty in the excitation energies is significantly smaller 
than the uncertainty in the total energies: apparently, 
the calculated total energies share a significant system- 
atic uncertainty, an overall binding uncertainty, which 
cancels when results are expressed as excitation energies. 
Furthermore, we see that both extrapolation methods 
agree very well with each other (within their error es- 
timates), and that the error estimates decrease as one 
increases the basis space. 

The extrapolation B leads to results for the two lowest 
excited states that are practically independent of the os- 
cillator parameter Ml, see Fig.[2j Also the bare and LSO 
renormalized NCSM results for these two states show 
very little dependence on Ml. These states are narrow 
resonances, and agree very well with experiment. 
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TABLE II: NCFC results for the 6 Li ground state E ga and excitation energies E x (in MeV) obtained in different basis spaces 
with JISP16. For extrapolations A and B we include in parentheses an estimate of the accuracy of the total energies; for 
the LSO(2) renormalized interaction, we p resent the spread in excitation energy for HQ variations from 12.5 to 22.5 MeV. 
Experimental data are taken from Ref. [la |. 





Extrap. A 


Extrap. B 


LSO(2) 


Extrap. A 


Extrap. B 


LSO(2) 


Experiment 


E{r,T) 


iVmax = 2 — 8 


iVmax = 4-8 


Nvxex = 6 


AT max = 10-16 


iVmax = 12-16 


iVmax = 14 


Energy 


Width 


£ 9S (l + ,0)! 


-30.9(0.6) 


-31.1(0.3) 




-31.47(0.09) 


-31.48(0.03) 




-31.994 


Stable 


■M3 + ,0) 


2.6(0.5) 


2.5(1.2) 


2.2-2.7 


2.56(0.04) 


2.55(0.07) 


2.53-2.55 


2.186 


24 • 10~ 3 


£*(0+ 1) 


3.6(0.6) 


3.5(1.2) 


3.3-3.7 


3.68(0.06) 


3.65(0.06) 


3.6-3.8 


3.563 


8.2 • 10~ 6 


E X (2+,Q) 


5.3(0.9) 


5.5(1.8) 


4.8-5.8 


4.5(0.1) 


4.5(0.2) 


4.8-5.0 


4.312 


1.30 


E x (2+,1) 


6.3(0.7) 


6.1(1.6) 


6.2-6.5 


5.9(0.1) 


5.9(0.1) 


6.0-6.4 


5.366 


0.54 


E X (1+,Q) 2 


6.1(1.7) 


6.6(0.3) 


7.1-8.5 


5.4(0.3) 


5.4(0.2) 


6.1-6.6 


5.65 


1.5 




extrapolated N = 12, 14, 16 

r max ' ' 
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FIG. 2: (Color online) NCSM results for the spectrum of 6 Li 
with LSO(2) renormalized (dashed) and bare (dotted) JISP16 
at iV m ax = 14, compared to NCFC extrapolations to infinite 
basis space (solid). Experimental data are from Ref. [l5l |. 



On the other hand, the three higher excited states have 
a much larger width, see Table fll] Our calculations for 
these broad resonances show a significant dependence on 
both Ml and iV max , in particular for the excited (1 + , 0)2 
state which has the largest width. The extrapolation 
B to infinite model space reduces but does not elimi- 
nate the Ml dependence. We further note that the Ml- 
dependence of these excitation energies is typical for wide 
resonances as observed in comparisons of NCSM results 
with inverse scattering analysis of a-nucleon scattering 
states [la ], and that the slope of the Ml dependence in- 
creases with the width of the resonance. This is consis- 
tent with the results presented in Fig. [5] the width of 
the (1 + , 0)2 state is larger than the width of the (2 + , 0)2 
state; the latter is larger than the (2 + ,l) state width. 
Thus, there appears to be a significant correlation be- 
tween the resonance width and the Ml dependence. The 
validity of the extrapolation to infinite model space is not 
entirely clear for these states. 

We noted earlier that the LSO renormalized interac- 



tion does not provide a monotonic approach to the infi- 
nite basis space for the binding energies and this prevents 
simple extrapolation. On the other hand, the excitation 
energies with the LSO renormalized interaction are of- 
ten quite stable with -/V max . However, it is important 
to realize that this does not necessarily mean that these 
excitation energies are numerically converged: they do 
depend on Ml. The dependence of the excitation en- 
ergies on Ml decreases slowly with increasing N max , as 
seen in Table ILT1 In fact, the excitation energies obtained 
with LSO(2) renormalized JISP16 are nearly the same as 
those obtained with the bare interaction, except at small 
values of Ml. as illustrated in Fig. [2j For most states, 
the NCFC provides better results for the excitation en- 
ergies, with less basis space dependence than the LSO (2) 
NCSM calculations in finite basis spaces. Nevertheless, 
we can employ the LSO procedure to obtain estimates of 
the binding and excitation energies in small basis spaces 
where there are no NCFC results. 

We summarize our results for the spectra of 14 F and 
14 B in Table IIIII The excitation energies are obtained as 
a difference between the extrapolated total energies of the 
excited state and that of the ground state (see Table UJ). 
The spectra are rather dense and the spacing between 
energy levels is smaller than the quoted numerical un- 
certainty, which is that of the extrapolated total energies 
of the excited states. However, as discussed above, we 
expect that for narrow resonances the actual numerical 
error in the excitation energy is (significantly) smaller 
than the error in the total energy. 

Figure [3J shows that different excited states can have 
very different convergence behavior. (Although we pre- 
sented in Fig. [3] the 14 B results, the behavior of the 14 F 
states is similar.) At N max = 8, there are five low-lying 
excited states; the excitation energy of these states de- 
pends only weakly on the basis space as iV max increases 
from 2 to 8. Then there are numerous higher excited 
states which depend strongly on the basis space: their ex- 
citation energies decrease rapidly with increasing iV max . 
Only after extrapolation to the infinite basis space do 
they appear at excitation energies comparable to the 
other low-lying excited states. We see a similar phe- 
nomenon in NCFC calculations of other nuclei. 
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TABLE III: NCFC results for the 14 F and 14 B excitation energies E x (in MeV). For extrapolations A and B we include in 
parentheses an estimate of the accuracy of the total energies; for the LSO(2) renormalized interaction, we present the spread 
in excitation energy for Ml variations from 12.5 to 22.5 MeV. Experimental data are taken from Ref. [171 ] - 







NCFC and NCSM ab initio calculations with JISP16 






Experiment 






14p 






14 B 








E(J n ,T) 


Extrap. A 


Extrap. B 


LSO(2), iV max = 6 


Extrap. A 


Extrap. B 


LSO(2), AT max = 6 


r 


Energy 


E x (l~,2) 1 


0.9(3.9) 


1.3(2.5) 


1.4-2.2 


1.1(3.5) 


1.4(2.8) 


1.4-2.3 


(1- 


) 0.654(0.009)" 


£„(3-,2)i 


1.9(3.3) 


1.5(4.6) 


1.0-1.8 


1.7(2.9) 


1.4(4.6) 


1.0-2.1 


(3- 


) 1.38(0.03) 


E x (2-,2) 2 


3.2(3.5) 


3.3(3.5) 


3.3-3.7 


3.3(3.1) 


3.3(3.8) 


3.5-3.8 


2" 


1.86(0.07) 




Q O/Q 0\ 


2.8(4.8) 


2.0-2.6 


3.1(2.9) 


2.7(4.8) 


9 D— "3 1 


IA~ 
I 4 


















? 


[2.32(0.04)]' 
















? 


2.97(0.04) 


E x (l-,2) 2 


5.9(3.5) 


5.4(4.6) 


5.8-6.4 


5.9(3.1) 


5.5(4.8) 


5.7-6.4 






E x (0-,2) 


5.1(5.4) 


5.8(1.0) 


5.8-10.5 


5.5(4.8) 


6.1(1.4) 


4.9-10.4 






E x (l~,2)3 


6.2(4.8) 


6.3(2.8) 


7.2-11.5 


6.4(4.3) 


6.4(3.1) 


6.1-11.3 






E x (2~,2) 3 


6.4(4.6) 


6.3(3.4) 


7.3-10.9 


6.9(4.1) 


6.7(3.6) 


6.6-10.9 






E x (3-,2) 2 


6.9(4.2) 


6.4(4.6) 


7.6-10.6 


7.0(3.7) 


6.5(4.7) 


6.4-10.5 






E x (5~,2) 


8.9(3.5) 


7.9(5.9) 


9.2-11.0 


8.8(3.1) 


7.8(5.9) 


8.5-10.8 







"Updated from Ref. O. 

'The existence of this state is uncertain (see Ref. [17|). 




FIG. 3: (Color online) Negative parity 14 B spectrum obtained 
with JISP16 at fixed Ml = 25 MeV in successive basis spaces, 
and extrapolated to infinite basis space using Extrapolation 
B. Experimental data are taken from Ref. [13 ] . 



The dependence on H£l varies considerably over the 
excited states as seen in Fig. SJ The lowest five excited 
states have a weak dependence on Ml, whereas the higher 
excited states depend strongly on it. We expect our re- 
sults for these higher excited states to have a larger nu- 
merical error than our results for the lower excited states 
with the weaker Ml dependence. Furthermore, in anal- 
ogy to the excited states in 6 Li discussed above, we expect 
these higher states to be broad resonances. Interestingly, 
the high-lying J" = 5~ state has a relatively weak Ml 
dependence (compared to states with similar excitation 
energy); it is also less dependent on A^ max , and may cor- 



FIG. 4: (Color online) NCSM results for the negative par- 
ity spectrum of 14 B with LSO(2) renormalized (dashed) and 
bare (dotted) JISP16 at -/V max = 6, compared to NCFC ex- 
trapolations to infinite basis space (solid) from N metx = 4—8. 
The most reliable HQ value for this extrapolation method is 
at Ml = 25 MeV for all states depicted. Experimental data 
are taken from Ref. [l3| . 



respond to a narrower resonance. 

Note, the conventional wisdom suggests leading con- 
figurations for the ground and 5 lowest-lying levels of 
Up ( 14 B) to be formed by the P3/2 neutron (proton) and 
S1/2 or CZ5/2 proton (neutron). Other low-lying states, 
with the exception of our low- lying 5 _ state, involve P1/2 
and/or d 3 / 2 single-particle states. 

In Fig. @] we can also see that the excitation energies 
obtained with LSO(2) renormalized JISP16 are nearly 
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the same as those obtained with the bare interaction, at 
least at N max = 6. Note that the NCFC results differ 
significantly from the bare and LSO(2) results, in par- 
ticular for the higher excited states with a strong iV max 
dependence; these extrapolated results also tend to have 
a somewhat weaker dependence on Ml than the results in 
finite basis spaces, and are expected to be more accurate. 

Some of the excited states in 1 B were observed experi- 
mentally. Unfortunately, the spin of most of these states 
is doubtful or unknown. Overall, the calculated exci- 
tation energies appear to be too large when compared 
with the experimental data; in particular our result for 
the excited 2~ state, the only excited state with a firm 
spin assignment, is about 1.5 MeV above the experimen- 
tal value. However, the spin of the lowest five states 
agrees with experiment, except for the 2~ and 4~ being 
interchanged, assuming that the tentative experimental 
spin assignments are correct. We do not see additional 
states between 2 and 3 MeV, but this could be related to 
the fact that all our excitation energies appear to be too 
large. Also, given the strong dependence on Y max of the 
higher excited states, it is not unlikely that these states 
will come down further with increasing basis space. 

We performed the first theoretical ab initio study of 
the exotic proton-rich nucleus 14 F which has yet to 
be observed experimentally. Using the NN interaction 
JISP16, we presented a prediction for the 14 F binding 
energy that is supported by comparing our NCFC re- 
sults with experimental data for the binding energies of 
the mirror nucleus 1 B and the neighboring nucleus 13 
obtained within the same approach in the same basis 
spaces. 

We extended our NCFC extrapolation techniques to 
evaluate excited states, and validated this method by ap- 
plying it to excited states in 6 Li. The obtained spectrum 
for 14 B agrees qualitatively with the limited data, and 
we made predictions for the spectrum of 14 F. More defi- 
nite information about the excited states in 14 B would be 
helpful. It would also be very interesting to compare our 



predictions for the F binding energy and spectrum with 
the experimental data that are anticipated soon. Sig- 
nificant differences between our predictions and the ex- 
perimental results would indicate deficiencies in the NN 
interaction, JISP16, and/or the role of neglected 3NF's. 
This would inform future research efforts and, with the 
inclusion of additional unstable nuclei in the analysis, 
aid in the eventual determination of the underlying shell 
structure evolution. 

Although NCSM calculations with LSO(2) renormal- 
ized interactions generally give reasonable results for 
the binding energies and spectra in small basis spaces, 
they do not improve systematically with increasing basis 
space. In particular for JISP16 we find that the results 
for the bare and the LSO(2) renormalized interaction ba- 
sically coincide for Y max > 8, both for total energies and 
for excitation energies. It would be worthwhile, although 
it is a major undertaking, to evaluate the effects of the 
induced three- and four-body interactions, which should 
improve the accuracy of the LSO renormalized calcula- 
tions. Without a thorough study of those effects however, 
we prefer the NCFC approach, based on extrapolations 
of NCSM results with the bare interaction, at least for 
JISP16. 
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